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The charged particle multiplicity distribution is one of the first measurements that ALICE 
will be able to perform. The knowledge of this basic property at a new energy is needed to 
configure Monte Carlo generators correctly with the aim of understanding the background 
of other, especially rare, processes including new physics. It allows to study the scaling 
behaviour and to verify model predictions. The unfolding of the measurement is a non-trivial 
task due to the finite precision and acceptance of the detector. Solutions are based on % 2 
minimization or iteratively using Bayes' theorem. Both approaches to unfold the spectrum 
are presented. Furthermore, the capabilities of the SPD fast OR trigger are shown that 
enable physics at very high multiplicities. 



I. INTRODUCTION 

A Large Ion Collider Experiment (ALICE) [l[ will take data in p+p (^/s = 14 TeV) and Pb+Pb 
(V s nn = 5.5 TeV) collisions. The charged particle multiplicity distribution P(N), i.e. the probabil- 
ity P for a collision with the charged particle multiplicity N, is one of the first measurements that 
can be performed. The knowledge of this basic property at a new energy is needed to tune Monte 
Carlo generators with the aim of understanding the background of other, especially rare, processes 
including new physics. It allows to study its scaling behavior and to verify model predictions. 

The multiplicity distribution of p+p(p) collisions at low energy is well-described by KNO scaling 
0|, that is broken at SPS (Super Proton Synchrotron) energies ( v / i = 540GeV) 0]. Other models 
and descriptions use negative binomial distribution (NBD) Q, a two component approach 0] or 
consider multiple parton interactions [6j. Only few predictions for LHC's energy realm exist, e.g. 
0] (based on the Quark-Gluon String Model) and 0] (based on the combination of two NBDs). 

ALICE can measure the multiplicity distribution in a broad region of phase space. The region of 
—3.4 < r) < 5.1 can be accessed combining the measurement of the Forward Multiplicity Detector 
(FMD) and the Silicon Pixel Detector (SPD). Figure [T] shows the coverage in 77 of the ALICE 
subdetectors. The analysis presented in these proceedings focuses on the measurement in the 
central region, using the information from the SPD in \r/\ < 1.4 and the measurement with the 
Time Projection Chamber (TPC) (\r]\ < 0.9). The latter can be used as a cross-check. 



II. ANALYSIS 



From all measured events only events that have a reconstructed primary vertex are considered 
for the analysis. This allows to select events from a given vertex range. Furthermore, the vertex is 
used to form tracklets in the SPD and to select good-quality TPC tracks. 

For each event the tracks that have a reasonable quality are counted. This step results in a 
measured multiplicity spectrum. Subsequently, this spectrum has to be unfolded in order to obtain 
the best estimate of the true multiplicity distribution (called true spectrum in the following) . This 
is non-trivial due to the fact that events with different true multiplicities contribute to the same 
measured multiplicity. The unfolding leads to the true spectrum - the multiplicity distribution of 
charged primary particles - for the events that have been triggered and have got a reconstructed 
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FIG. 1: The Figure shows the acceptance in 
rj of all ALICE subdetectors overlayed with a 
dN c h/dri prediction by Pythia 9] for y/s = 
14TeV. All subdetectors have full coverage in az- 
imuth except the ones marked with an asterisk. 



FIG. 2: The Figure shows the detector response of 
the SPD, i.e. reconstructed tracklets vs. the num- 
ber of generated primary particles in |ry| < 2.0. 
Only events with a z vertex position in ±10 cm 
are considered. 



vertex. Finally, this spectrum needs to be corrected for the bias introduced by the vertex recon- 
struction requirement and the trigger. 



III. UNFOLDING 



The response of the detector can be described by a matrix R mt . It gives the conditional prob- 
ability that a collision with a true multiplicity t is measured as an event with the multiplicity m. 
The response matrix is created from the full detector simulation (see Figure [5] for the response 
of the SPD). 

Given a true spectrum T, the measured spectrum M can be calculated by M = RT. Therefore, 
given a measured spectrum, the true spectrum is formally calculated as follows T = R~ l M. 
However, R^ 1 cannot be calculated in all cases, because R may be singular. Even if R can be 
inverted, the result obtained contains strong oscillations (due to statistical fluctuation caused by 
the limited number of events used to create the response matrix) [j^ . Alternatives to simple 
inversion are Bayesian unfolding and \ 2 minimization, introduced in the following paragraphs. 

Bayesian Unfolding is an iterative procedure based on Bayes' theorem which describes the 
definite relationship between the following two conditional probabilities: the probability of an event 
A (m measured particles) conditional on another event B (true multiplicity t) and the probability 
of B conditional on A. Applying the theorem on the problem at hand results in [13, [l3|: 

R tm = JfZ tP l (1) U t = -J2M m R tm (2) 

V 

R is the response matrix and P the a-priori distribution of the true spectrum which can be 
chosen to be a flat distribution or like in the present analysis the measured spectrum. defines 
i?, which is in literature also referred to as smearing matrix. Subsequently, {2]| with R and the 
measured spectrum M allows to obtain U. tt is the efficiency to detect an event with a given true 
multiplicity. It is different from unity due to the vertex reconstruction and trigger efficiency. In the 
next iteration, U is used as a-priori probability P. Optionally a smoothing is applied that reduces 
the influence of high-frequency fluctuations (e.g. a sliding average over three multiplicity bins). 
After about 10 - 100 iterations the result converges, U approximates the true distribution. The 
left panel of Figure [3] shows the result of a Bayesian unfolding after 100 iterations with smoothing. 




The strength of the smoothing as well as the number of iterations are free parameters that have 
been evaluated in detail in (T3 |. 

% 2 Unfolding is an approach to unfold the measured multiplicity distribution by the mini- 
mization of a x 2 function. Using the response matrix, this function gives a measure of how well a 
(guessed) unfolded spectrum describes the measured spectrum. A minimization program is used to 
find the unfolded spectrum that has the least \ 2 . With e m denoting the error on the measurement 
of M and U being the guessed spectrum, a suitable function is: 

P(U) is called regularization or penalty term. It only depends on the guess of the unfolded 
spectrum (and not on R or M). It is a constraint that favors a certain shape of the unfolded 
spectrum and thus prevents the previously mentioned fluctuations. f3 determines the weight that 
is given to the regularization with respect to the first term that governs the agreement with the 
measurement. A higher (3 increases the value of the first term in ^ because U is more constrained. 
Its optimal value needs to be evaluated (see [3]); a reasonable value of (3 adjusts the two terms 
in © in a way that the introduced bias is negligible compared to the statistical error of the 
measurement. 

Many possibilities exist for the choice of the regularization P: these range from "just" requiring 
a smooth function to preferring a certain shape of the distribution. Generally, not too much knowl- 
edge about the distribution should be implied, otherwise the result is likely to look very similar to 
the required function and may mask details of the measurement. Various regularization functions 
have been evaluated in [l4j , a suitable one is to require the least curvature: 

^)-E(f)' = E( P "-^ + [fw )' ,4) 

An unfolding result is shown in the right panel of Figure [3] that uses regularization (j4]) with 
(3 = 10 4 . 



IV. HIGH MULTIPLICITY TRIGGER 



The multiplicity spectrum can be measured up to very high multiplicities by exploiting the SPD 
fast OR trigger. This trigger is based on the number of fired chips in the two layers of the SPD (radii 




of 3.9cm and 7.6cm from the beam line). A chip is fired if at least one particle deposits energy in 
one of the chip's attached pixels. In total 1200 chips produce a trigger at a rate of 10 MHz. Trigger 
thresholds on the number of fired chips select high-multiplicity events; several thresholds can be 
configured simultaneously. The left panel of Figure 0] shows for the first layer < 2) the number 
of fired chips vs. the event multiplicity together with a trigger treshold, set at 150 chips; the right 
panel shows a multiplicity distribution obtained by simulating a combination of the minimum bias 
trigger and two high multiplicity triggers. It reaches up to 400 charged particles, thus enabling to 
study e.g. the abundance of hyperons like the S at high multiplicities which provides insight in the 
hyperon scaling behavior with multiplicity. 



V. SUMMARY 



Two unfolding methods have been presented to obtain the multiplicity distribution. Either of 
them reproduces simulated spectra successfully. For the measurements, both methods should be 
used, as consistency check. Optimizations in both methods have been evaluated and the results 
described in these proceedings. More details and the study of the systematic effects that arise 
during this measurement can be found in A study performed on ALICE'S high multiplicity 
trigger shows that this trigger promises interesting insight into the physics at high multiplicities. 
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